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Abstract 

In a well-known paper [Ij V.C. de Andrade, L. C. T. Guillen and 
J.G. Pereira defined a conserved gauge current hj^^^', however they stated 
that: "TLis is, we believe the farthest one can go in the direction of a 
tensorial deEnition for the energy and momentum of the gravitational 
field. The lack of local Lorentz covariance can be considered as the 
teleparallel manifestation of the pseudotensor character of the gravi- 
tational energy-momentum density in general relativity . . " . Well, we 
believe that they stopped just less than an inch before giving such a ten- 
sorial definition, and furthermore that the resulting energy-momentum 
tensor has zero trace and can be made symmetric, as a matter of fact it 
is just j^^^^ , and together with the energy-momentum tensor of material 
fields it obeys a natural conservation equation for teleparallel manifolds. 
Some important consequences are obtained specially in the last section 
concerning the possibility of explaining the acceleration of the universe 
expansion without any need of a cosmological constant. 



1 The energy— momentum tensor 

We follow most of the conventions and notations of references PP and 
However one difference is that we stick to consider the lagragian density of 

o 

the gravitational field to be proportional to (as it is usually accepted) 

instead of being proportional to —R as is done in [Ij . Any further difference 
will be made it explicit. Anyway let us at least remember that greek indices 
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are used for the coordinate holonomic quantities and latin indices are used for 
non-holonomic ones. 

One way of viewing the kind of Riemannian spaces in which teleparallel 
theories are formulated is by taking as departure point the hypothesis that 
physics somehow establishes a canonical, smooth, path-independent isomor- 
phism between the tangent spaces of any two points of the manifold and hence 
that we may take some orthonormal, but otherwise arbitrary reference ba- 
sis, which we will call and refer all vectors at all points to that basis (of 
course, wc identify Ua with its preimage at any given point). Parallel transport 
can then be introduced as meaning to transport keeping constant components 
with respect to this basis, then Cartan covariant derivative is just the varia- 
tion with respect to this reference basis expressed, for example, in terms of the 
coordinate basis. Mathematically, this has the consequence of accepting only 
parallelizable manifolds as physically meaningful, this is nothing more than 
just a topological condition on the sort of Riemannian spaces we deal with. In 
particular the coordinate vectors d^{x), no matter which sets of coordinates x 
we take, should be expressible in terms of the reference basis: 

= h%{x)Ua (1) 

Let us also remember that the relationship between the Lcvi-Civita covariant 

o 

derivative V due to the symmetric riemannian connection and the "Cartan" 
covariant derivative V of the Weitzenbock connection is given by the difference 
between the Christoffel symbols of both covariant derivatives: 

tiv — HI/ + ij^u (2) 
Where pn is the contorsion tensor given by: 

Now, as a first point in the discussion we must clarify what can be considered 
a conservation equation for the energy-momentum tensor within teleparallel 
theory. Suppose we had some (symmetric) energy-momentum tensor and 
let be the components of some vector which is Cartan covariant constant. 
Such vectors do exist because any linear constant combination of vectors of 
the reference basis Ua is Cartan covariant constant. S^^w^^ represents the flow 
of the component of energy-momentum in the direction of the four-vector w, 

o 

so \7j^(^S^'^w^) = expresses the conservation of such w component. However 
we can write this as: 

= w''(V.V-V^%) (4) 
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If we want conservation of energy-momentum in all directions, then it must 
hold: 

V.S/ - S;T\, = (5) 

o 

This is not the same as the condition VuS^^ = 0. If we substitute the Cartan 
covariant derivative by its classical counterpart, we reach another expression: 

V.^/ - S/K%, = (6) 

Of course, these equations reduce to the zero divergence condition in Minkowski 
spaces, however the important point is that (for second rank tensors in telepar- 
allel spaces) they represent the correct generalization of the zero divergence 
condition of Minkowski space. 

o 

As it is well known Einstein's tensor G^'^ verifies Vi/C^'^ = 0, so from 
Einstein's equation we do not get a conservation equation for the energy- 
momentum tensor of matter alone. For the energy-momentum tensor T of 
matter, accepting it to be a symmetric tensor, we rather get: 

= V.r^ ' =^ (V.T^ - "T^p) + % ^T\, = (7) 

Comparing it with ([5]) we see it is not quite exactly the same, a further term 
is present. Let us remember that we are dealing with a lagrangian density 
A for the gravitational field which in teleparallel theories might in general be 
written as: 

A = Kg{aiAi + asAa + asAs) (8) 

where: 

Ai = g'^T^xT^^ A2 = g'-'T'^^T^^^ (9) 
As = r*r,,„ = ^ (10) 

The "default" values: ai = 1, 02 = —1/2, 03 = —1/4, produce a lagrangian 
density equivalent to general relativity, meaning that for those values the dif- 

o 

ference between R and (ciiAi + a2A2 + ctsAs) is just a total divergence. As a 

„ o , o 

matter of fact, adding 2g°''^\/aT fsx one obtains R. However, for the moment, 
we are not going to fix the value of those coefficients, we want to keep open 
the possibility of choosing other values for them. 

One way of obtaining the gravitational energy-momentum tensor is to di- 
rectly reproduce the reasoning which in classical mechanics leads to the con- 
servation of energy and write (we follow the standard work P], which as a 
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matter of fact is just the first step of applying Noether's method considering 
the translational invariance of the lagrangian, see also [S], section 2): 



d{Ah) _ d{Ah) dh\ ^ d{Ah) dh\^^ 



The field equations can be written as: 



1 

h' 



'd{Ah) d 
dh"„ 



h- 



dA 



r 



And taking them into account one inmediately is led to: 



1 d 



hdx'f 



h h' 



dA 



dh'' 







(11) 



(12) 



(13) 



One would like to identify the term within the round brackets in last equation 
with the energy-momentum tensor, but it is not a tensor, so the idea is to 
decompose it into tensorial and non-tensorial terms, hence we write previous 
equation as: 

1 d 



hdx^ 



HQ,' + n; 







(14) 



Where is the tensorial part and A^^'' is a non-tensorial term. Expressing 
the partial derivative of the Q tensor as a Cartan covariant derivative one 
arrives inmediately to the following equation: 



dN„^ 
dx"' 



(15) 



Now the problem is to eliminate the non-tensorial terms of this equation. Of 
course, it must be possible to do so, because you cannot have an equalty be- 
tween entities which transform in different ways (you may put the first two 
terms at one side and the other terms at the other side). The needed calcu- 
lations for eliminating, or transforming, those non-tensorial terms are a bit 
cumbersome, but in the appendix it is indicated how it can be shown that by 
taking Q^^ as: 



+ a2 
+ as 



Kg {a, [2 ((?^^T^,TV-^7"^2^ 
2 ((7'^^T%,T^„^-(7'^"TVr''. 



a rp'y 



4^7^^^7""^7p.T^^T^J - [ai Ai + asAs + a.A^] 5;} (16) 



one arrives to the conclusion that 



dN„^ 



per /T- u 



(17) 
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The expression for might seem strange at first sight, but it happens to be 
exactly -j^^': 

^ = 7 ^ = -Q ^ (18) 

Hence we are led to the result: 



V,j7 - 3;T%, + % -"T^, = (19) 

So taking into account both equations (fT9l) and ([7]) one gets: 

V,(T^ ' + j/) - (T/ + j7)T%. = (20) 

Which has exactly the form of equation ([5]) and so it can be interpreted 
as just expressing the conservation of total energy-momentum: the energy- 
momentum of the material fields 7^ plus the energy-momentum of the grav- 
itational field j^''. It also clarifies the meaning of the term T^^T" This 
term specifies the energy-momentum interchange between the gravitational 
field and the material one. It is the term which prevents energy-momentum 
of the gravitational field or energy-momentum of the material field from being 
conserved separately by themselves. The interpretation of as the correct 
energy-momentum tensor for the gravitational field can be further underlined 
if one rewrites the field equations as: 

and compares it with the equations for the electromagnetic field in Minkowski 
space written as: 

d f dAe\ r 



(22) 



Informally, it is usually accepted that this last equation says that currents are 
the sources of the electromagnetic field. Then, in the same sense, the previous 
equation might be interpreted as saying that energy-momentum of both the 
material fields and the gravitational field is the source of the gravitational field. 

However, this is not the final word about the tensor we need, because 
j^^ has the uncomfortable characteristic that in general it is not symmetric. 
If we accept the default values for the coefficients 01,02,03 then there is an 
easy way out of the problem: Einstein's equations are a total of ten equations, 
however they do not completely determine the metric tensor. Diffeomorphisms 
comprise the gauge freedom in general relativity (see [6] pg. 438): any two 
solutions which are related by a diffeomorphism represent the same physical 
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solution. Now, if instead of considering the metric tensor as the final solution 
of a gravitational problem, we ask a complete solution of such a problem to be 
given by the specification of the sixteen functions h°'^ which give the coordinate 
basis vectors in terms of the arbitrary constant reference basis, then we have 
some further freedom, because once the metric tensor is given, we may have 
several "square roots" h"'^ which produce the same metric tensor: There are 
sixteen arbitrary h""^ functions and only ten independent conditions imposed 
by the metric tensor. We have lots of "gauge freedom", so let us use part of 
that freedom to decree that the antisymmetric part of the energy-momentum 
tensor jj^ should be zero. This "gauge condition" amounts just to a set of six 
equations, because in a four dimensional space an antisymmetric second rank 
tensor has only six independent components. So, although it is a crude way of 
counting degrees of freedom, we have increased by six the number of unknown 
functions when substituting the metric tensor as solution by the h^^, but we 
have also added six additional equations, so we expect not to have changed 
the "gauge freedom" . The gauge condition can be written as: 

= = [g^^T'^p, - g'-T'fs,) g^^T^ax - \ {g^'^T- - g'^'T^^,) (23) 

Of course it is a covariant condition: true in one coordinate system means 
true in all, so we are not limiting the set of coordinates in which the theory 
is formulated: we are not imposing conditions on the sort of diffeomorphisms 
which might be used. As a matter of fact, it has previously been argued 
that teleparallel theories may have too much gauge freedom (see [7], [8]) and 
that they suffer from a problem of non-predictability of torsion. Although a 
formal proof should be investigated, we clearly expect this gauge condition 
to fix such problems. At least, the introduction of this condition invalidates 
the reasoning supporting such assertions, because clearly these additional six 
equations have not been taken into account when studying the predictability 
of torsion. Furthermore: being an algebraic condition on the torsion tensor, 
not every boundary condition is acceptable, because the boundary condition 
must also obey the gauge condition. 

Accepting such a gauge condition, the energy-momentum of the gravita- 
tional field turns out to be symmetric, which just means that j"''^ = And 
furthermore it is inmediate to check that it has zero trace. Needless to say it is 
a perfectly covariant local definition of energy-momentum for the gravitational 
field. 

There is one further point which merits some comment. Teleparallel theo- 
ries have some degree of freedom in the way the coefficients are chosen. How- 
ever it is only for the case in which we obtain the teleparallel equivalent to 
general relativity when the resulting equations are symmetric (we obtain Ein- 
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stein's tensor). So it is only in this case in which we have the freedom to 
impose that the energy-momentum tensor of the gravitational field should be 
symmetric (as there is no other possibility in empty space). So this condition 
eliminates the rest of possibilities for the coefficients. 

We may even generalize: if we restrict teleparallel theories to just this case, 
as it has been said before, we obtain symmetric equations (Einstein's tensor), 
so what would happen if the energy-momentum tensor of a material field were 
not symmetric?. The most natural answer in this hypothetical case would be 
to change the gauge condition so that the antisymmetric part of the energy- 
momentum tensor of the gravitational field just cancels the antisymmetric part 
of the energy-momentum tensor of the material fields, and only the symmetric 
part plays a role in Einstein's equation. 

2 The energy content of homogeneous, isotropic 
universes. 

The energy content of homogeneous, isotropic universes has already been com- 
puted in other papers (see [3j, for example). However, although our results are 
quite similar, there are some points which should be noted: previous works 
have used pseudo-tensors and of course have never taken into account the 
gauge condition. Usually they need to integrate over a space section. Given 
that we are dealing with a isotropic, homogeneous universe, if we were really 
working with a true gravitational energy density, one would expect the energy 
density to be homogeneous, and if the total energy turns out to be zero, it 
is quite unintuitive that it is not zero at every point. Of course, one answer 
is that usually authors are not working with the true energy density of the 
gravitational field. 

To calculate the energy content of the three cases of null, positive and neg- 
ative curvature, we must first find "square roots" of the respective metrics. 
Not every square root is acceptable, we must also impose the gauge condition: 
the energy-momentum tensor derived from them should be symmetric. How- 
ever for these problems we must not only obtain symmetrical tensors, but also 
diagonal, otherwise we would have some preferred direction in space. It must 
be noted that most calculations in this section have been done with Maple© 
9.5 running on Ubuntu Linux. 

The easiest case is the one of zero curvature. We use a "cartesian" coordi- 
nate system with coordinates ct, x, y, z, and we postulate the following matrix 
of gravitational potential vectors (the role played by the coordinate vectors is 
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similar to that of vector potentials): 



diag(l,a(i),a(i),a(t)) 



(24) 



Using this potentials, the metric is just the very well known diagonal metric 
of fiat space g^^ = diag(l, —a^{t), —a'^{t), —a^{t)). The energy-impulse tensor 
for such a space is: 



Where the dot signals ordinary differentiation with respect to ct. The first 
point which deserves attention is that energy density is negative, so it seems 
that there is at least a known field whose energy density takes negative values. 
It is also purelly "kinetical" in this case: it is proportional to the square of 
the speed at which a{t) changes, and the minus sign tells us that absortion of 
(positive) energy will decrease this speed. 

Let us consider first the case of a dust- filled universe. We know that for 
such a case a{t) = Cdt'^^^, so the gravitational energy density is proportional 
to — which, when multiplied by oc to take into account the increase in 
volume, just gives constant energy: dust does not contribute to any variation 
of energy of the gravitational field, it does not interchange energy with the 
gravitational field. 

Consider now the case of a universe filled with just radiation, being p its 
energy density. It is well-known that in such a case the solution for a(t) is of 
the form a{t) — Crt^^^. So the gravitational energy density is also proportional 
to —t~^ (it is the square of a logarithmic derivative, so no matter the exponent 
it will be proportional to —t~^), which when multiplied by a^{t), to take into 
account the increase of volume, gives the result that energy of gravitational 
field changes as — which is an increase and which is just the rate needed to 
compensate the rate at which energy of radiation decreases: pa^ is constant, so 
pa^ decreases as a"^ oc The absortion of positive energy from radiation 

just decreases the rate at which universe expands. In the dust-filled case, the 
decrease in speed is just to compensate the increase in volume, so that total 
energy is the same. As energy of light is absorved by the gravitational field, its 
"kinetic" energy increases (it decreases its "speed"). A radiation dominated 
universe expands at a slower rate (i^^^) than a dust filled one (i^/^): absortion 
of energy decreases its speed. 

Even more clear, for the flat imiverse Einstein's equation can be used to 
calculate the energy-momentum tensor of matter: 




(25) 




(26) 
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Looking at the first term we see that the energy density of the matter fields is 
just the same (but positive) as the energy density of the gravitational field, so 
that total energy density is zero. This is a fact which also happens in the case 
of positive curvature. 

For positive curvature universe the way to get an acceptable square root of 
the metric is to remember that is parallelizable, so it is easy to get three 
ortonormal vectors which are tangent to it and then take them as the spatial 
part of the reference basis at each point. Of course, by an acceptable square 
root wc mean one which renders a symmetric gravitational energy-momentum 
tensor. A diagonal square root, as it is usually taken, docs not. The previous 
idea can be summed up by saying that we postulate an /i"^ matrix given by: 
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a(i)c(e) -a{t)s{^)c{^)s{e) -a(t)s'^(^)s'^{e) 

a(t)s(e)c($) a(t)s(*) + c(<I')c(e)c(*)) a(t)s(*)s(e) (s(*)c(e)c($) - c(*)s($)) 

a(t)s(e)s(*) a(t)s(*) (c(*)c(e)s(*) - s(*)c(*)) a(t)s(*)s(e) (c(*)c(*) + s(*)c(e)s(*)) 

(27) 

Where s = sin and c = cos. This matrix leads to the spherical metric: 

g,^ = diag(l, -a'{t), -a\t)s\^), -a\t)s\^)s\Q)) (28) 

The energy-momentum tensor of matter is given by: 

+ 1 ^1 ^ „ 2da + d^ + 1 , , 

= 6k,— ^ ^i' = ^ = = 2k, (29) 

CL (Jj 

The gravitational energy-momentum tensor is given by: 

^o = -K— ^ jI = J2 = jI = ^^9—^ (30) 

(Jj (Jj 

So we get into the same situation in which the total energy is zero. 

It is a bit more difficult to find an acceptable square root for the negative 

curvature case. As a matter of fact, in the acceptable solution which has been 
found, the potential vectors cannot be so neatly separated into spatial and 
temporal parts, however they do lead to the correct metric and to a diagonal 
energy-momentum tensor. The solution found for matrix h"'^ is: 



/ ch(^) a(i).s/i(*) 

s/i(*)c(e) a{f})chi;^)c{e) -a{t)sh{'^)sie) 

s/i,(*)s(e)c($) a{f^ch{^)s{e)c{<^>) a{t)sh{^')c{e)c{(^) -a(t)sh{^).s{e).s{<S>) 

\ s/i(*)s(e)s($) ait)ch{^)s{Q)si<S>) a(t)s/i(«')c(e)s($) a(f)s/i(«')s(e)c(4>) / 



(31) 
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Where ch = cosh and sh = sinh. It is inmediate to check that the hyper- 
bohc metric is obtained with these potentials: 

g^^ = diag(l, -a^{t), -a'^{t) smh^{^), -a^{t) sinh^(^) sin2(e)) (32) 

The energy-momentum tensor of matter is given by: 

tO a ~ ^ Ti ^2 tS o ^^'^ + a? — I 
— 6/«c/ ; 7^ = ^2 = = ^Kg — ; (33) 

While the gravitational energy-momentum tensor is given by: 

fo = -6-.^^ jI = jI = jI = 2S^^^ (34) 

(Jj LL 

But we do not reach the same conclusion: the total energy is not zero. This 
is somewhat unexpected: if the total energy were zero in the three models, 
it will agree quite well with the idea of originating from a state of zero total 
energy. However this case is different and so, unless some other interpretation 
is found, it raises the question of whether, although mathematically possible, 
it is physically reasonable: why should universe begin in a state of energy 
different from zero?. On the other hand, a flat universe has the minimum 
matter density compatible with zero total energy. 



3 Last comments. 

We have already argued that equation expresses the conservation of energy- 
impulse, so let us write that equation as: 

V = - S^^T\, = (35) 

We have also seen in equation ([7]) that Einstein's equation implies: 

O^r^ ^ = -T^ ^T%^ (36) 

Where T is the energy-momentum tensor of matter fields. Let us suppose we 
are dealing with a perfect fluid in an isotropic homogeneous universe. Let us 
consider for example the case of fiat universe. The energy-momentum tensor 
can be written in such case as: 

%y = diag (^p, pa^ , pa^ , pa^^ (37) 

Where p = p{t) is the mass-energy density, p = pit) the pressure and the a? it) 
factors come from the metric. The right hand of equation (136|) can be easily 
computed and turns out to be: 

(^-3%(t),0,0,o) (38) 
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The first thing which stands out is that if p{t) ^ then in an expanding 
universe, mass-energy of the material field (by itself) is not conserved. We 
have seen such a behaviour when we considered the radiation-filled universe 
before: a positive pressure means that the gravitational field drains the positive 
energy from the "material" field. As a matter of fact we have also seen the 
case p = in the dust-filled universe and there was no energy interchange. 
Let us turn to the other possible case. Suppose there is some "spontaneous 
matter emission" process, then if matter is created from the gravitational field, 
pressure must be negative. Of course, we do not know what exactly to put in 
the left hand, and the matter-emission process must have very low probability 
of ocurrence because otherwise it would have already been detected. Let us 
just put a small "constant" A in the left hand side of equation flHUjl . mainly 
because we have no better guess, then we may write: 

3p(t) = -A^ (39) 

So although the process may have very low probability of ocurrence, the neg- 
ative pressure increases with the expansion of the universe. It may of course 
overcome the mass term in the equation which determines the acceleration of 
the expansion of the universe (see for example [6J pg. 97): 

3- = --^[p(t) + 3p(t)] (40) 
a AKg 

From that moment, positive acceleration sets in. We do not need to have a 
cosmological constant to explain it. In fact we cannot consider this A as a 
constant, it may depend on the strength of the gravitational field, and also 
even in case we accept the possibility of matter-emission processes, their rate 
must compensate the energy absortion rate of the gravitational field from elec- 
tromagnetic radiation. We would need a quantum theory of gravitation to be 
able to calculate X{t). However, we may get an idea of the order of magnitude 
of A by considering zero the acceleration, taking the Hubble constant Hq = a/a 
to be 70 (km/s)/Mpc, and taking the density of the universe to be the critical 
one 2 X 10~^^kg/m^. We get X ^ 5 x 10"^''kg/(m^s). This is the order of 
magnitude of the rate at which matter is created at the expense of the gravita- 
tional field (of course, it says nothing about what sort of particles are created). 
There is nothing to prevent gravitational field from falling even further down 
in energy levels (as it is the usual objection to negative energies), only that the 
rate is extremely slow. Supposedly, a quantum theory of gravitation should 
be able to explain this rate. Anyway it sets an experimental test for any such 
a theory: see if within lowest order perturbation you can obtain something 
similar. 
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A The calculation of Q. 



First place we will need the explicit form of Einstein equations written in terms 
of the torsion tensor and metric tensors. This has been done quite a number of 
times, and we only put explicitly and directly in terms of these tensor. After 
working out the calculations implicit in equations (fT^ . Einstein's equations 
can be written as: 



Kg 



ai 



2g"'V^T"xa + 5:(2Ao - Ai) + [2V^T^^ - 2V^T^\ 



+ 



-ATP^-'T^p^ - AT\pT„ + 2T'^-"T.„, - Ag^^V^T^^"" + Ag^^l (41) 



To accomplish the elimination of nontensorial elements in equation ffT^ let us 
do the needed calculations for each of the three parts of the lagrangian density. 
For Ai = g^'^T'^axT'^p^, one has that: 



And one may take: 

Kn 



-2 (g'^T\xT%, 



(42) 



2 {g^^T\^T%, - T\^T^'^) - A^6^ 



and 



Kn 



2 (g^^T' 



a A J- 1^1(3 



P7 ' 



(43) 



(44) 



The subindex 1 of, for example, lA^^ just refers to the fact that the term comes 
from Ai, it is not any sort of spatial index. So after some work, which may 
involve using the condition that the curvature of the Weitzenbock connection 
is zero for simplifying some expressions, one finds that: 



dx"/ 



- {2g^''V^T\a + 5l{2Ao - Ai)) F'^ 

For A2 = g^'^T'^pxT'^af, one gets: 
dA2 



-2 (g'^'T^xT^^u - g'^'T'^BxT^.u) - A,6Z 



(45) 



(46) 



i/,7 
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And then one may take: 



and 



-2 ((7"^T"^Ar^M"-^?''^Vr'^ 



(47) 



(48) 



So after some second work, one finds that: 



Finally for A3 = g^'^g'^'^ gp^T^ axT" one has: 
(9A 



And then one may take: 



Kn 



^g^^g'^^gpaTPo^xT^p.-A^dl 



and 



Ag^^g'^^gp^TP^xT' 



pl3 



(49) 



(50) 



(51) 



(52) 



So after some third work, one finds that: 



1 



-4F%, (-T^^^^g, - T\pT/P + ^T'^^^T.^A - V^T,'^'^ + Ag^^) (53) 

So taking into account Einstein's equation fHTl) . one arrives to the conclusion 
that: 



d(N„^) 



(54) 
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